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ABSTRACT 

Baryon acoustic oscillation (BAO) is a powerful probe on the expansion of the universe, shedding light 
on elusive dark energy and gravity at cosmological scales. BAO measurements through biased trac- 
ers of the underlying matter density field, as most proposals do, can reach high statistical accuracy. 
However, possible scale dependence in bias may induce non-negligible systematical errors, especially 
for the most ambitious spectroscopic surveys proposed. We show that precision spectroscopic redshift 
information available in these surveys allows for self calibration of the galaxy bias and its stochasticity, 
as function of scale and redshift. Through the effect of redshift distortion, one can simultaneously 
measure the real space power spectra of galaxies, galaxy-velocity and velocity, respectively. At rel- 
evant scales of BAO, galaxy velocity faithfully traces that of the underlying matter. This valuable 
feature enables a rather model independent way to measure the galaxy bias and its stochasticity by 
comparing the three power spectra. For the square kilometer array (SKA), this self calibration is 
statistically accurate to correct for 1% level shift in BAO peak positions induced by bias scale depen- 
dence. Furthermore, we find that SKA is able to detect BAO in the velocity power spectrum, opening 
a new window for BAO cosmology. 

Subject headings: cosmology: distance scale-large-scale structure of universe-theory 



1. INTRODUCTION 

The universe is the largest laboratory for fundamental 
physics. The zeroth order descriptions of the universe 
include the overall expansion rate H(z) and its integral 
form, the distance as a function of redshift z. Type la 
supernovae (SNe la), as promising cosmological standard 
candles, allow the measurements of cosmological distance 
and have enabled th e discovery of the late time accelera- 
tion of the universe (R iess et al. 1998: Pcrlm utter et alj 
1999). This discovery has profound implication for fun- 
damental physics, leading to either a dominant dark en- 
ergy component with equation of state w = P '/ p < —1/3 
or significant deviations from general relativity around 
Hubble scale. Ongoing and planned supernova surveys 
will be able to significantly narrow the range of w and 
hopefully clarify the role of the cosmological constant in 
our universe. Meanwhile, other independent probes are 
indispensable, to reduce statistical errors, break param- 
eter degeneracies and to cross check for possible system- 
atical errors. The baryon acoustic oscillation (BAO) is 
such a probe dBlake fc Glazebrook 2003: Hu & H aimanl 
120031 : ISeo fc EisensteiiJl2003L 120051 fBlake et al.ll2006| ). 

Prior to the epoch of recombination, baryons and pho- 
tons are tightly coupled. Fluctuations in the baryon- 
photon fluid propagate as acoustic waves at speed com- 
parable to the speed of light. These acoustic waves im- 
print in the cosmic microwave background, a snapshot 
of the photon fluid at the epoch of decoupling. They 
also imprint in the baryon fluid through coupling be- 
tween photons and baryons and later in the overall mat- 
ter density field, through gravitational coupling between 
baryons and dark matter particles (|Hu fc Sugivamal 
U996tlEisenstein fc Hulll998l) . These BAO features have 
been detected and measured to high precision in CMB 
and have put strong constraints on the geometry of the 



universe (|Spergel et al.ll20"07"[) . In the large scale struc- 
ture (LSS) of the universe, these features are signifi- 
cantly suppressed by dark matter to a level much weaker 
than that in CMB. Nonetheless, they are detectable, as 
shown in measurements of correlation fu nction and power 
spectrum of 2dF and SPSS galaxies (ICole et all 120051: 
lEisenstein et al.ll2005t IPercival et al.ll2007t ). The comov- 
ing scales of BAOs can be predicted accurately from well 
understood linear physics. Thus we are able to convert 
the measured angular separations of these features on 
the sky to the cosmological angular diameter distance, 
free of many astrophysical uncertainties. Furthermore, 
since we can measure their separations in redshift, we 
are able to directly measure H(z), which is more sen- 
sitive to w than the distance. For these merits, BAO 
is widely accepted as a robust and powerful cosmologi- 
cal stand ard ruler to reyeal th e expansion history of the 
universe (jAlbrecht et al.ll2006l ). 

All existing proposals on BAO measurements target 
at density field of various tracers of the LSS. Most are 
galaxy spectrosc opic redshift surveys, i ncluding LAM- 
OST 1 WieeleZ (iGlazebrook et al J 12001. BOSS 2 , WF- 
MOS (|Bassettet a ll 120051), HET DEX 3 , ADEPT 4 and 
the s quare kilometer array (SKA, lAbdalla fc Rawlingsl 
2005). Besides, there are proposals to measure 
BAO in the 21cm background, both at l ow redshifts 
(IChang et all 120071) and a t high redshifts (|Mao fc Wul 
I2007t IWvithe et al.l I2007D, i n the Lyman-alpha for- 
est dMcDonald fc Eisensteinl 1 2 0071 ) and in SNe la 
(|Zhan et al.l 120081 ). Sub-percent level statistical accu- 
racy can be reached in the most ambitious surveys, 

1 http://www.lamost.org/en/ 

2 http://cosmology.lbl.gov/BOSS/ 

3 http:/ /www. as. utexas.edu/hetdex/ 

4 http: / /www7. nationalacademies.org/ssb/BE_Nov_2006_bennett.pdf 



2 



such as SKA, classified as a stage IV projec t by the 
dark energy task force (|Albrecht et al.l 120061 ). How- 
ever, possible sources of systematical errors could pro- 
hibit us to exploit full advantage of these measurements. 
Some, for example those induced by gra vi tational lensing 
(IVallinotto et al.ll2007tlHui et al.ll2007al|b1: lLoVerde et al 



20071) are correctable (|Zhang fc CheniT2007l : iHui et al 



2007b). Nonlinear evolution in the matter density field 



is also likely correct able, either by improvement i n the- 
oretical calculation (ICrocce fc Scoccimarrol 120061 120071: 
iMatarrese fc Pietronil l2007allbl; iMcDonaldl \200% or in 
data analysis ( Eisenstein et al.M2007b[ K However, the 
clustering bias of these tracers with respect to the under- 
lying matter density field, is more elusive to handle. The 
scale dependence of the bias changes the overall shape of 
the power spectrum and thus shifts the BAO positions. A 
good thing is that, to the first order accuracy, BAOs, as 
periodic sharp features on top of the smooth power spec- 
trum, are robust against bias scale dependence, which 
is likely more smooth. For this reason and others, sys- 
tematical errors caused by possible scale dependent bias 
are likely sub-dominant, comparing to the statistical er- 
rors, for BAO measurements of existing surveys and pro- 
posed medium cost surveys. On the other hand, the 
most ambitious surveys such as SKA will have about two 
orders of magnitude increase in survey volume and one 
order of magnitude increase in redshift over SDSS. For 
these surveys, statistical errors in the distan ce and H(z) 
meas urement will be reduced to below 1% (jBlake et al.1 
2006). At this stage, shifts in BAO peaks induced by 
even rather mild scale dependence in galaxy bias may 
become non-negligible. For discussions on th e effec t 
of bias to BAO, refer t o e.g. ISeo fc Eisenstein! <l2005f): 
lEisenstein et al.l (l2007ah: iGuzik et all (l20q7[): iHuff et all 
(|2007fk ISmith et al.l (|2007afk lAngulo et all (|2008[ )~ 

Fortunately, spectroscopic redshift information, which 
is required to measure H{z) along the line of sight 
through BAO, allows for self calibration of galaxy bias 
through the effect of redshift distortion. Spectroscopic 
redshift surveys are able to measure Pg (k) , the 3D power 
spectrum of galaxies in redshift space, to high precision. 
Galaxy peculiar velocities make Pg anisotropic in the 
3D wavevector k space. This redshift distortion effect 
is often treated as a source of error. However, it can 
be rendered into valuable source of signal instead, given 
sufficiently accurate understanding of redshift distortion. 

Through the anisotropy in redshift space, one can re- 
construct the real space velocity power spectrum Pg, the 
velocity-galaxy cross correlation power spectrum P g g and 
the galaxy power spectrum P g in a rather model indepen- 
dent way at cosmological distances. Here 8 is the nor- 
malized divergence of peculiar velocity. Since galaxies 
have rather low mass, they are sub-dominant in gravita- 
tional field at scales > lOMpc. At scales of BAO (much 
larger than 10 Mpc, even for the third peak), gravity is 
the only force accelerating galaxies. For these reasons, 
galaxies can be well approximated as test particles, fol- 
lowing the motion of dark matter particles. Thus the 
large scale galaxy velocity field is a faithful tracer of the 
overall matter velocity field. Pg and its relation with the 
matter power spectrum P m can then be predicted from 
first principles, free of uncertainties in understanding of 
galaxy formation. Comparing P g , P g g and Pg, we are 



able to measure galaxy bias and its stochasticity. This 
allows for self calibrati on of galaxy bias in BAO cosmol- 
ogy In a pioneer work, iPenl (|1998f ) proposed to measure 
galaxy bias and its stochasticity through the quadrupole 
and octupole of Pg . This work extends his by a full power 
spectrum reconstruction, for which we are able to utilize 
information contained in all moments. 

On the other hand, BAOs also exist in the velocity 
field. Pg reconstruction thus opens a new window of 
BAO measurement. This approach is free of problems 
of galaxy bias, although it may not be able to reach the 
same statistical accuracy as that of the density field. 

This technique can be incorpo rated into existing BAO 
analy sis met hods, such as that of[Seo fc Eisenstein! (1200 
120071) and iBlake fc Glazebroo kl (12003D; iBlake etal. 
(|2006l) ; iPercival et alj ( 20071 ). as straightforward post 
processing. We will explain this issue later. 

SKA is able to detect billions of galaxies over half sky 
to redshift z > 3 through neutral hydrogen 21cm emis- 
sion of these galaxies and measure precise spectroscopic 
redshifts without extra cost. For this reason, we choose 
SKA as our primary target for investigation. We present 
the reconstruction procedure in demonstrate the fea- 
sibility of SKA to self calibrate galaxy bias in show 
that BAOs in the galaxy velocity field are detectable 
through SKA in Sj4] and discuss in $5] 

2. VELOCITY RECONSTRUCTION THROUGH REDSHIFT 
DISTORTION 

Peculiar motions of galaxies imprint unique signa- 
tures in the redshift space galaxy power spectrum Pf (k) , 
which takes a general form (e.g. IScoccimarrol (|2004| )). 



P 9 S (k) = [Pg(k) + 2u 2 P g g(k) + U 4 Pg(k)] F 



H 2 (z) 



1) 



Here, P g , P g g, Pg are the real space power spectra of 
galaxies, galaxy-velocity and velocity. is the comoving 
peculiar velocity divergence divided by —H(z). a v is the 
ID velocity dispersion; and F(x) is a smoothing function, 
normalized to unity at x = 0, determined by the veloc- 
ity probability distribution, u — k\\/k is the cosine of 
the angle of the k vector with respect to radial direction; 
This unique directional dependence has enabled success- 
ful simultaneous reconstruction of the three power spec- 
tra fr om 2dF and SDSS galaxies (|Tegmark et all 120021 
12004 . We a dopt a simple minimum variance estimator, 
developed in I Zhang et al.l |2007), to quantify accuracies 
of this approach. 

For each kj in the given k bin, we have a measure- 
ment of Pg, which we denote as Pi. The unbiased mini- 
mum variance estimator of the band power of Pi a \ (a — 



9, 
A 



is P = Y.W ( t a) P l , where W ( t a) 



F, 



(a) 



+ Xf'uf) and Fi = F(k Ul a v /H). a i =P l + l/n g 
is the rms fluctuation of P, , where n g is the galaxy num- 
ber density. The Lagrange multipliers Aj°2 3 are given 
by = (1,0,0) • A -1 , \ {ge) = (0,1/2, 0) • A^ 1 and 
A( e ) = (0, 0, 1) • A -1 . A(") are orthogonal to each other. 
The 3x3 matrix A is given by 



E 



2(m+n-2) F» 

2a 



i m, n = 1, 2, 3 



(2) 
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The corresponding error in each power spectrum is 

One can check that it has the right scaling that a p oc 
a/yN, where N is the number of independent modes 
and a is some average of a%. This relation can be further 

simplified and we have <j [ p 9) = \f \[ g) /2, a p 9d) = y/>^J~4 
and a p 9) = 

Combining these power spectrum measurements, we 
are able to measure the galaxy bias and its stochasticity 
through 

9 \Pe) \Pm) \heory V 

Here, (Pg / P m ) theory is predicted from theory. In linear 
regime, this value is f~ 2 , where / = (d\nD/d\na) 2 and 
D is the linear density growth factor, r is the cross cor- 
relation coefficient, w hose deviation from unity is a m ea- 
sure of stochasticity ()Penlll998l : iDekel fc Lahavlll999t ). 

This reconstruction requires the angular diameter dis- 
tance Da and H (z) as input to convert the observed an- 
gular and radial separation into k. For numerical results 
presented in this paper, we just take Da and H (z) of the 
fiducial cosmology as input. For real data, a convenient 
and self consistent procedure to carry out this reconstruc- 
tion is as follows. First, one can measure the distance 
and H (z) by the usual met hods (|Blake fc GlazebroolJ 
120031: ISeo fc Eisensteinl l2003ft . The measured distance 
and H(z) are then used as input for the reconstruc- 
tion of the power spectra, which are in turn applied to 
check for the consistency of scale independent bias. This 
can be done iteratively. Alternatively, one can allow for 
scale dependen t bias and its stochastic i ty in BAO anal- 
ysis method of ISeo fc Eisensteinl (|2003l I2007D . Through 
simultaneous multi-parameter fitting, galaxy bias (am- 
plitude, scale dependence and stochasticity) is automat- 
ically self calibrated. 

Errors in Da and H(z) propagate into the power spec- 
trum reconstruction, through two effects: a constant 
fractional shift in fc and a wrong determination in the 
k direction u. Obviously, the overall shift in fc does not 
affect the measurement of b g and r. However, a wrong 
determination in u does. P g g and Pg are determined by 
the derivatives of P g (u) with respect to u. We can show 
that the induced errors in the power spectra have only 
weak dependence on k through the slope of the power 
spectra. Thus the main effect of errors in Da and H(z) 
is to bias the overall amplitude of b g and r, instead of 
introducing new scale dependence. For these reasons, it 
suffices to neglect this kind of error source. 

An implicit simplification adopted in this reconstruc- 
tion is that a v is given. In reality, at BAO relevant scales, 
the function F(kua v /H) ~ 1. Thus our simplification is 
justified here. For real data, to improve the accuracy of 
reconstruction, one should treat free parameter 

to be marginalized over. We do not expect any major ef- 
fect on the forecasted errors in the reconstructed power 
spectra. 

Hereafter we will adopt the linear theory to carry out 
numerical calculations. For real data, nonlinearity must 



be taken into account in BAO analysis. For example, 
nonlinear evolutio n in density field and velocity field 
are no t the same (|Bernardeau et all l2002t ISmith et al.1 
2007b). As a result, (Pg / P g ) theory is no longer equal to 
f~ 2 and is no longer scale independent. If unaccounted, 
we may obtain a false scale dependence in b g . How- 
ever, this should not be a real problem, due to three 
reasons. (1) The reconstruction does not rely on assump- 
tions of linearity. (2) The bias is defined with respect to 
the real (thus nonlinear) power spectra. (3) Nonlinear 
evolution in the velocity field and matter field can be 
calculated from first principles. As long as we replace 
the linear theory version f~ 2 with the actual (nonlinear) 
(Pg / P m ) theory in real data analysis, nonlinearity does not 
bias the result. 

3. SELF CALIBRATION OF GALAXY BIAS 

In this section, we choose SKA as the primary tar- 
get for investigation. For the redshift range 1 < z < 2 of 
most interest, the SKA survey speed (A e ff/T sys ) 2 xFOV= 
2 x 10 10 m 2 K~ 2 deg 2 . 5 With this specification, SKA is 
able to detect several billions of 21cm emitting galaxies 
in 10 year survey. Throughout this paper, we adopt the 
galaxy surface number density as 20 galaxies per square 
arc-minute, corresponding to 2.4 billion over three quar- 
ter of sky. The exact number is hard to estimate, due 
to poor knowledge on these galaxies at high rcdshifts. A 
good thing is that the exact number is not necessary for 
the estimation presented in this paper, due to the dom- 
inance of cosmic variance over shot noise, especially at 
scales larger than the second baryon acoustic peaks. 

To proceed, we adopt the fiducial cosmology as a flat 
ACDM cosmology, with Cl m = 0.268, Q, b = 0.044, f_ A = 
1 - J l m , h = 0.72 and cr 8 = 0.77 (jSpergel et all 
120071). We adopt the tra nsfer function fitting formula 
of lEisenstein fc Hul (|1998l ). We stick to the linear per- 
turbation theory and neglect the nonlinear evolution in 
the matter density field, since it is predictable and cor- 
rectable. The velocity power spectrum is fixed through 
Pg = f 2 P m - Here, / = d\nD/d\na, D is the linear den- 
sity growth factor and P m is the matter power spectrum. 

To demonstrate the feasibility of self calibration, we 
introduce scale dependence and stochasticity in galaxy 
distribution. As a reminder, the galaxy bias b g (k,z) is 
defined by b 2 (k,z) = P g (k, z)j P m (k, z). The stochas- 
ticity of galaxy distribution is characterized by the 
cross correlation coefficient r(fc, z), defined through r 2 = 
P 2 g/(P g Pg)- Deterministic bias will have r = 1. This is 
likely true at very large scales. Proceeding to smaller 
scales we may expect r < I, because galaxy forma- 
tion depends on not only the local density, but also the 
environment. We adopt simple forms of parametriza- 
tion: b g = b g o + db g /dk\k=ok = b g o(l — fc/fcb) and 
r = 1 + dr I dk\k=ok = 1 — k/k c . Here, b g o is the galaxy 
bias when k — > 0. We further adopt fcb = 4/i/Mpc and 
k r = 10/i/Mpc. The actual scale dependence of bias 
can be much more complicated, with terms oc fc, oc fc 2 , 
oc 1 /P m (k), or even mo re complicated terms (refer to, 
e.g. ISmith et al.1 l2007at ) . We choose the above simple 

5 http: / /www. skatelescope.org/pages/Preliminary_spccifications_v2. 4.pdf. 
This survey sp eed is much higher t h an wh a t adopted in preivou s 
estimations of lAbdalla fc Rawlingsl p005h : IZhang fc Pent (|20Q6h . 
resulting in much higher sensitivity to detect high redshift galaxies. 
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forms of parametrization and values for several reasons. 
(1) These forms of parametrization are just the Tay- 
lor expansion of b g and r at k = to first order. De- 
spite their simplicity, they are rather general to describe 
small deviations from scale independence. (2) For the 
values adopted, the galaxy bias decreases by less than 
5% over the first three BAO peaks. This scale depen- 
dence has the same sign as, but smalle r amplitude than, 
what predicted for low mass galaxies dSheth fc Tormenl 
fl99l iGuzik et al.ll2007t iSmith et all l2007aft . Thus the 
assumed scale dependence likely represents a conserva- 
tive lower limit of this source of systematical errors. (3) 
Such level of scale dependence in b g may cause non- 
negligible systematical errors. It shifts the positions of 
the second and the third peak by ~ 1%. This level of sys- 
tematic shift in BAO positions could overwhelm the sta- 
tistical uncertainties in SKA BAO measurements. 6 For 
these reasons, the assumed bias scale dependence and 
stochasticity are suitable to demonstrate the feasibility 
of the self calibration technique. 

We remind the readers that, these parametrizations 
of b g (k,z) and r(k,z) only serve as input of the fidu- 
cial model. The reconstruction of P g , P g g and Pg does 
not rely on any assumptions on b g and r. Parametriza- 
tions on b g and r can certainly improve over the statisti- 
cal errors in the reconstruction. However, inappropriate 
parametrizations could bias the reconstruction. For this 
reason, we do not explore such approach in this paper. 

For SKA, P g can be reconstructed to high precision 
(Fig. [T]) and the associated fractional statistical er- 
rors are very close to that of corresponding P*. Ac- 
curacies in P g g reconstruction are about a factor of 5 
poorer, but still impressive. A further factor of ~ 3 
degradation occurs in Pg reconstruction. This is what 
we expected. The average contribution of Pg to 
is Pg(u 4 ) = Pg/5 ~ P m il 1 m 2 (z)/5 and that of P g g is 
2P g g(u 2 ) ~ P m 2n° m 6 (z)/3. The realistic galaxy bi as is 
likely b g > 1/2 (Fig. [2] and further reading in Uingl 
EMI). We thus have Pg(u A ) < 2P g g(u 2 ) < P g . Fur- 
thermore, we rely on the directional dependence of Pg to 
extract sub-dominant components P g g and Pg, resulting 

in 4 e) /P 6 > O^lPgg > 4 9) /Pg. 

Nonetheless, the reconstructed power spectra (espe- 
cially Pg) have the statistical accuracy to detect 5% or 
even smaller variation in the galaxy bias over the first 
three peaks, with large statistical significance (Fig. [T]). 
This 5% scale dependence causes 1% shift in BAO peak 
positions. Thus the self calibration technique is able to 
reduce systematical error in BAO peak positions induced 
by scale dependent bias to smaller than 1%. 

However, despite its promising capability, the self cal- 

6 The exact level of systematical errors induced by scale de- 
pendent bias varies wit h BAO analysis methods. The method of 
Sco & Eiscnstcin (2003) performs the global fitting to the galaxy 
power spectrum and utilizes the cosmological information in the 
overall sh ape of the power spectrum . On the other hand , the 
method of IBlake fc GlazebrooH |(2003T ): IPercival et al.l p007l ) dis- 
cards such information. It fits the overall shape of P g with a 
smooth reference function P re f. It then fits P g /P re { by a decay- 
ing sinusoidal functiion, in which the BAO information is encoded. 
P rc f can capture and thus filter away some of, if not most of, the 
bias scale dependence. Thus this method is less affected by the 
scale dependence in galaxy bias. But the residual could still bias 
the distance measurement. 
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k (h/Mpc) 

Fig. 1. — Forecasted accuracies in reconstructed power spectrum 
Pgi Pg9 ar >d Pg- The estimation is based on SKA ten year survey 
over three quarters of sky and galaxies at 1.5 < z < 2.5. We 
plot the ratio of three power spectra with respect to Pg (Q(, — * 0) 
and normalized the ratios to be equal at k — > 0. Error bars are 
forecasted adopting b g o = 1. Pg is much harder to measure than 
P g and P g g, so the self calibration is fundamentally limited by 
the associated error in Pg. Nonetheless, SKA is able to test and 
correct for several percent level scale dependence in bias across 
the first three peaks. Furthermore, with such accuracy, BAOs are 
detectable in Pg. Accuracy in Pg is sensitive to b g , as shown in 
Fig. E] 



ibration procedure is fundamentally limited by accuracy 
in our understanding of redshift distortion. Systematic 
deviation in the modelled Pg from the actual one can 
propagate into the estimation of b g and thus can bias 
the self calibration of galaxy bias. A good news is that, 
to measure the scale dependence and perform self cali- 
bration, we are not necessarily confined to BAO scales. 
As seen from Fig. [TJ we are able to reconstruct the three 
power spectra to k ~ 0.3/i/Mpc and even smaller scales. 
k = 0.3/i/Mpc is still in the linear regime at z = 2. This 
means that we have a factor of ~ 2 larger (and still safe) 
dynamical range to measure the scale dependence in b g . 
Furthermore, variation in b g over this larger dynamical 
range is likely larger. Both alleviate the accuracy require- 
ment for modeling the redshift distortion (e.g. Eq. [T]). 
For the simple model we adopted, this means to mea- 
sure 8% variation in b g over k < 0.3/i/Mpc instead of 
measure ~ 5% variation over k < 0.2/i/Mpc. Still, the 
modeling of the redshift space power spectrum must be 
accurate to much better than 2 x 8% = 16% percent in 
order to render the associated systematical error negligi- 
ble. This requirement is challenging, but not impossible, 
given the rapid development in simulations and analyti- 
cal calculations. We refer the readers to works on mod- 
eling of redshift distortion, through perturbation theory 
(IScoccimarrcj|2004lMatsubaral|2007l ). halo model (| White) 
120011: lSeliakH2001t lTinkerll2007f) and N-body simulations 
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Fig. 2. — The SKA galaxy selection threshold (top panel) and 
the corresponding bias (bottom panel). SKA galaxies are selected 
by the galaxy neutral hydrogen mass Mm, while the bias is deter- 
mined by the total galaxy mass M. Lack of robust measurements 
on /hi = Mni/M, we plot the galaxy bias for several possible 
values of /hi- 
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Fig. 3. — Baryon wiggles in Pg of 0.5 < z < 1.5 galaxies. The 
accuracy of Pg reconstruction is sensitive to b g . Smaller b g causes 
stronger directional dependence in P g and weaker statistical fluc- 
tuations in P g and P g g and thus improves Pg reconstruction. SKA 
is able to detect numerous galaxies with weak clustering at z ~ 1. 
Pg reconstructed through those galaxies can be sufficiently accu- 
rate for BAO measurement in the velocity field. 



(IKang et alJl200l iTinker et al.l l2006). 

4. DETECTING BAO IN THE VELOCITY POWER 
SPECTRUM 

BAOs in the density field inevitably show up in the 
gravitational potential field, as directly seen from the 
Poisson equation. This featured gravitational field accel- 
erates galaxies and thus imprints BAOs in the velocity 
power spectrum Pg. At scales of BAOs, the galaxy veloc- 
ity filed should be unbiased with respect to the velocity 
field of the underlying matter, since these galaxies are 
test particles responding only to gravity at these scales. 
For this reason, we are able to calculate Pg from first 
principles. Since in principle there is no ambiguity in 
precision calculation of Pg , BAO cosmology based on Pg 
is free of the issue of galaxy bias. It is thus worthy of 
detecting BAOs in Pg. 

In reality, BAO measurements through reconstructed 
Pg from redshift distortion are challenging, since accu- 
racies in Pg are much worse than that in Pg or P g , as 
explained in last section. For b g o = 1, the fractional error 
in Pg is a factor of ~ 10 higher than that in P g (Fig. [1]). 
It thus requires ~ 100 times increase in survey volume 
in order to reach the same statistical accuracy obtained 
from the density field. For this reason, BAO can not be 
detected in Pg by those low redshift spectroscopic sur- 
veys such as LAMOST and BOSS. But SKA can (Fig. Q] 
fc[51), so as the Hubb le sphere hydrogen survey (HSHS) 
((Peterson et al.ll2006ft and ADEPT. For SKA, the statis- 
tical accuracy in Pg is about 3% around the second and 
the third peaks (for bin size Afc = 0.01/i/Mpc), slightly 
larger than that of the P g measurement in the combined 
sample of SDSS main galaxies and luminous red galaxies 



(|Percival e t al. 2007). SKA extends much larger redshift 
range than SDSS. Higher redshift means more k modes 
suitable for BAO detection (less severe smear of BAO fea- 
tures by mode-coupling in nonlinear regime) and higher 
sensitivity to dark matter and dark energy. Thus BAO 
cosmology through the recovered Pg may gain over cur- 
rent SDSS BAO measurement through P g . 

Accuracy in Pg reconstruction is sensitive to the galaxy 
bias b g . A smaller bias causes a stronger directional de- 
pendence and enables better Pg reconstruction. At the 
first three peaks of BAO, shot noise is in general sub- 
dominant for SKA. This allows us the luxury to select 
those low mass, weakly clustered galaxies for better BAO 
measurement in the velocity field. A practical question 
is that, how strong are SKA galaxies clustered. Obser- 
vationally, we know that galaxies in local universe de- 
tected through 21c m line emission are weakly clustered 
(|Mever et al J 12007). However, we are lack of observa- 
tions on b g of neutral hydrogen mass selected galaxies 
other than in the local universe. For this reason, we just 
present a rough estimation here. 

SKA galaxies are selected by their neutral hydrogen 
m ass Mhi- We update the calculation of lower Mm limit 
in I Zhang fc Pen! ((2006) with new SKA specifications. In 
a ten year survey, Mhi at 3<r selection threshold is shown 
in FigO To calculate the associated bias of these galax- 
ies, we need to convert Mhi to the galaxy total mass M. 
We define the neutral hydrogen fraction /hi = Mm/M. 
The value of /hi is poorly known. But we are able to de- 
rive its lower limit. The tot al neutral hydrogen density 
is H i ^ 3.5 x 10 -4 today (|Zwaan et al.l l2005f) and in- 
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creas es by a factor of ~ 5 toward z ~ 2 — 3 (|Peroux et al.1 
2003). Since most baryons and dark matter are not in 
galaxies while most neutral hydrogen atoms are in galax- 
ies, we expect that /hi 3> f^m/^m- For several choices 
of /hi, w e plot the corr esponding b a in Fig. [2l We adopt 
the bg-M relation from lJingl (|19981 ). whose fitting formu- 
lae improves over previous analytical results, especially 
for low mass halos. Many of SKA detected galaxies at 
z = 1 are likely well below 10 11 M Q . These galaxies are 
indeed weakly clustered, with bias as low as b g ~ 0.6. 
These galaxies are very likely numerous, if the shape of 
neutral hydrogen mass function does not change at low 
mass end, from local universe to z ~ 1. Selecting only 
these galaxies then does not result in significant increase 
in shot noise. What we gain is a factor of ~ 2 improve- 
ment in Pg reconstruction and comparable improvement 
in BAO measurement, comparing to the case of b g o = 1 
(Fig. EJ). 

5. DISCUSSIONS AND SUMMARY 

Up to now we only discussed the feasibility to perform 
self calibration of galaxy bias and to measure BAO in the 
velocity power spectrum reconstructed from redshift dis- 
tortion of galaxies. This technique can be applied to any 
tracers of the large scale structure with precision redshift 
information. This includes the pseudo 21cm background 
composed of unreso lved 21cm emitting galaxies at z ~ 1 
(|Chang et al.ll2007l) and the 21cm background at higher 
redshifts (|Mao fc Wdl2007t iWvithe et al.ll2007lh 

It is also feasible to measure BAOs in the velocity field 
reconstructed from redshift distortion. A fundamental 
obstacle must be overcome is the statistical fluctuations 
associated with the reconstruction technique. Since sta- 
tistical fluctuations in P g and P g g propagate into Pg re- 
construction, statistical errors in Pg reconstructed are a 



factor of ~ 10 times larger than its cosmic variance (Fig. 
[J). To reach the cosmic variance limit of Pg and improve 
significantly on BAO measurements in velocity field, new 
techniques of velocity measurement are needed. Kinetic 
Sunyaev Zel'dovich effect of galaxy clusters is a possibil- 
ity. It is also possible to do precision peculiar velocity 
measurement through milli ons of SNe la with s pectro- 
scopic redshifts at z ~ 0.5 (iZhang k Chenll200l . BAO 
measurement brings extra scientific benefit for such sur- 
veys. 

Pg reconstructed through the galaxy spectroscopic sur- 
veys contain valuable information of cosmology, besides 
the information of the expansion of the universe con- 
tained in the BAO. Combining the BAO measurements 
through galaxy spectroscopic surveys and Pg, we are able 
to infer the initial fluctuations of the universe, the struc- 
ture growth of the universe and the geometry of the uni- 
verse simultaneously. This potentially very promising 
and powerful application will be addressed in a compan- 
ion paper. Pg and P g g, when combined with gravitational 
lcnsing, also have important applications in probing the 
nature of gravity at cosmological scales ()Zhang et al.1 
2007, 2008). Finally we emphasize that all these appli- 
cations rely on precision modeling of redshift distortion. 
Tremendous work is required to meet the requirement of 
precision cosmology. 
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